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Considering the earthquake susceptibility of structures located in seismic 
regions, the question arises as to what level oj protective measures should 
be provided in order to achieve a certain degree of reliability against possible 
damage. To address this question, engineering risk and optimal design 
of structures located in a seismic area are studied. The basic concept is 
to obtain a tradeoff between the cost of providing a protective measure 
and the expected cost of earthquake damage. 

A simple mathematical approach is presented to determine the optimal 
earthquake intensity which the structure is designed to withstand. The 
objective is to minimize the total construction cost of the structure plus 
the expected cost of earthquake damage throughout the entire service life 
of the structure. For the case of deterministic structural resistance, and 
for structural response processes having Poisson {independent) crossings, 
an objective function is derived in terms of the building and earthquake 
variables. The optimal design intensity can then be determined by mini- 
mizing the objective function with respect to the intensity variable. The 
resulting equations are relatively simple and can be easily handled for 
numerical studies and sensitivity analysis. Generalizations of the results 
for nondeterministic structural resistance and for structural response 
processes different from those having Poisson crossings are also indicated. 

As an illustration of the proposed approach, a hypothetical building 
with realistic seismicity and structural parameters is analyzed for its 
optimal design earthquake intensity. The construction, damage and total 
costs are obtained in terms of the intensity variable. The implications and 
sensitivity of the results are also discussed. 

I. INTRODUCTION 

Structures constructed in seismic regions are required to function 
properly in a forcing environment characterized by random earthquake 
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occurrences and intensities. The seismic environment including the 
expected earthquake-magnitude levels and the corresponding frequency 
of occurrence for different seismic-risk zones was described previously. 1 
The study was based on a statistical analysis of nationwide seismic 
data and may be used as a guide for the development of seismic design 
requirements on a global basis. Under localized situations, however, 
the seismic requirement for structures that are expected to adequately 
withstand the earthquake environments should be based on cost- 
reliability studies. During an earthquake of given intensity, there exists 
a probability that the response of the structure is greater than its 
resistance capability and, therefore, a probability of damage to the 
structure. The cost associated with this probable damage may be 
referred to as the "earthquake risk cost". Increasing the design intensity 
of the structure reduces the probability of damage, but at an increased 
cost of construction. Therefore, an optimal design earthquake intensity 
can be determined by achieving an appropriate balance between the 
construction cost and the earthquake risk cost. 

This paper presents a new analytical approach to the determination 
of the economically optimal earthquake intensity or other design 
variables for structures. The construction and earthquake risk costs 
are expressed in terms of design intensity and other parameters re- 
flecting the earthquake and structural characteristics. Minimization 
of the total expected cost of the structure yields the optimum structural 
design intensity in terms of such parameters as estimated cost of 
earthquake damage, unit construction cost, expected earthquake 
duration, and statistics obtained from seismological data for the par- 
ticular site. 

II. ANALYSIS OF DESIGN INTENSITY MODEL 

2.1 Objective Function 

Consider a certain seismic region in which a structure is located. 
Let K c (i ) and K d (i ) represent the construction and the earthquake 
risk or damage costs of the structure respectively, both being functions 
of the design intensity i„ measured from I to XII on the Modified 
Mercalli scale. The function K c (i ) may be regarded as a monotonically 
increasing function of i„ , while the function K d (i ), as would be expected, 
is a monotonically decreasing function of i . 

The optimum design intensity i*, may be obtained as a trade-off 
between these two functions by minimizing the total cost 

K(i„) = K C (Q + K d {i ) (1) 
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Fig. 1 — Sketch showing cost of structure against design intensity of earthquake. 

as shown in Fig. 1. The function K(i ) in eq. (1) is the objective function 
and i„ is the decision variable. 

The construction cost, K c (i„), can be written as 



KM = A,(f + C) 



(2) 



where A f is the floor area of building, / is the building cost per unit 
floor area, and C is the cost of earthquake protection per unit floor 
area. The earthquake protection cost obviously increases with the 
protection level, such as i a in this case. Therefore, C = C(i„) is a mono- 
tonically increasing function of i„ . It should be noted that the deter- 
mination of the function C(i„) depends on the type of structure, method 
of design, and is greatly influenced by the designer's personal judgment 
and experience. A reasonable first approximation can be made assuming 
C(i„) is linear with a coefficient c. Under this assumption eq. (2) can be 
written as 



K e (i ) = A,(f + ri„) 



(3) 



where A f j equals a fixed initial cost of building, and c is the earthquake 
resistance cost of the building per unit floor area per unit, intensity. 
To determine the function K H {i u ), analytical procedures can be 
effectively employed utilizing some basic knowledge about the structure 
and the random forcing environment. Let N(t) be a random variable 
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representing the total number of earthquakes to occur in time t. Further- 
more, for k — 1, 2, • • • , N(t), let p k be the probability that the structure 
fails given, the /cth earthquake occurs, and let e k , an identically dis- 
tributed random variable, be the associated loss. Therefore, the total 
cost of earthquake damage, Z, is 

Z = £ e kVk . (4) 

Taking expectation (denoted by overbar) of both sides of eq. (4), one 
obtains 



Z = N(t)ep (5) 



in which N(t) equals the expected number of earthquakes (of all in- 
tensities) in time t, e equals the expected value of the random earth- 
quake loss, and p equals the mean failure probability of structure 
given an earthquake occurs. The present worth of the expected value 
of Z is the earthquake risk cost K d (i ): 

K d (i ) = Zg(t) = N(t)epg(t) (6) 

where g(t) is a discount factor. 

The quantities on the right-hand side of eq. (6) will be discussed next 
in terms of the related design parameters which fall into two categories: 
the earthquake parameters and the building parameters. The earthquake 
parameters include the regional seismicity constants, the earthquake 
magnitude (m), intensity (i), duration (l ), amplitude (a) of the ground 
motions, and the statistics of these quantities. The building parameters 
include the mass (p), stiffness (k), natural frequency (o>„), damping (£„), 
height (h), the resistance or strength (/) of the structure, etc. 

2.2 The Average Number of Earthquakes, N(t) 

The quantity A'(t) depends on and can be estimated from the regional 
seismicity. Earthquakes can be considered to be a series of events 
randomly distributed on a real line (representing time), and the sequence 
of original times \t„} forms a point process. 2 It is further assumed that 
the joint statistics of the respective number of shocks in any set of 
intervals are invariant under a translation of these intervals; this 
implies that {t n } is a stationary point process. The stationary point 
process generalizes certain aspects of renewal processes; in particular, 
the interval lengths r k = t k — t k - { between successive events need not 
be independently or identically distributed. 
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The simplest stationary point process is the Poisson process. In- 
tuitively, the process \t„) can be approximated as a Poisson process 
if it represents rare events. More rigorously, it requires that t k be 
independently and identically distributed and follow a negative expo- 
nential function. The main deficiency of the simple Poisson model is 
its inability to describe the aftershocks which are often triggered by 
a large main shock. However, for most practical engineering purposes, 
the simple Poisson model for earthquakes appears to be adequate. 
In practice, an engineer is concerned with the earthquake risk of struc- 
tures located in some specific geographic areas. The risk depends heavily 
on the statistics of large earthquakes in these areas, and the omission 
of small earthquakes or aftershock processes is relatively unimportant 
in terms of earthquake risk. 

If \N(t); t ^ 0} is assumed to be Poisson with a constant rate, a, 
then 

prob[AT(0 = n] = ~fe~ at (7) 

and 

N(t) = at. (8) 

The parameter a per unit time t can be determined from regional 
seismicity data. 1 

2.3 The Mean Failure Probability of Structure, p 

The quantity p depends on both earthquake and building parameters. 
The earthquake intensity "i" to which the structure is designed will 
be the only decision variable considered in this formulation and all 
other parameters are assumed known. Let Y(t) = {max | y(t) | ; t t 
[0, t„] \ , where, /„ equals the duration of the structural vibration which 
is assumed approximately equal to the duration of the earthquake 
ground motion, and y(t) equals the response parameter (displacement, 
velocity, acceleration, stress, etc.) of the structure. For an earthquake 
with intensity i, failure could occur when the resulting structural 
response y(t) equals or exceeds the actual resistance, x, of the structure. 
The corresponding failure probability, p(i), can be expressed by 

p(i) = prob [Y ^ .r | earthquake with intensity i has occurred]. (9) 

The quantity p(i) is a function of the random variable i representing 
the earthquake intensity whose probability density function j,(i) can 
be found in terms of the regional seismicity and earthquake source 



2214 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1972 

geometry. Since i, in the Modified Mercalli intensity scale, takes only 
discrete integer values from one to twelve, the mean failure probability of 
structure, p is given by 

v = E pV)Mi), (io) 

I - I o 

in which i u equals the design earthquake intensity. From eqs. (9) and 
(10) it is clear that earthquake parameters enter the formulation of the 
problem through p(i) and fi(i), while building parameters enter the 
formulation through p(i) only. 

The density function j,(i) can be derived from an expression obtained 
by Cornell* 1 for the distribution function of earthquake intensity i: 

F t (i) = 1 - y IV exp (-£*)• (11) 

The governing assumptions for cq. (11) are: 

(i) The earthquake magnitude m is a random variable with in- 
dependent and negative-exponential distribution function 

F M (m) = 1 - e~ 0m , m ^ m (12) 

where m is some magnitude small enough, say 4, that events of lesser 
magnitude may be ignored by engineers, and /3 is a constant the inverse 
and inverse square of which represent the mean and variance of earth- 
quake magnitude m > respectively; 

(it) The intensity attenuation law is given by 

i = c, + c 2 m — c ;i In r (13) 

in which c, , c 2 and c : , are regional seismicity constants, 1 and r, the 
focal distance in miles, is the random variable representing the distance 
from the structural site to the location of an earthquake source on the 
fault line. 

(Hi) The earthquake has a line source (fault line) of length I (in 
miles) with uniform distribution. 
The parameters T and .7 are given by: 

r = exp [^(l + in.)] (14) 

J = 2 /"" - y$ „ , y = P C f-l (15) 



+ For example, c, / = 1, 2, 3 are semi-empirical constants on the order of 8, 1.5, 
and 2.5, respectively, for firm ground in Southern California. 4 
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in which d = min r and r — max r. 
It follows from eq. (11) that 

dFj(i) 0TJ ( j8 .\ Mlft 

m- *--£-■» rsv' (16) 

2.4 Determination of the Failure Probability of Structure, p(i) 

To find the structure's failure probability p(i), as denned by eq. (9), 
it is necessary to specify the failure mechanisms of the structure. It 
is also necessary to establish stochastic models for the response pa- 
rameters y(t) and the corresponding resistance x of the structure. For a 
linear and deterministic structure which is assumed to experience no 
plastic deformations and the properties of which are governed by given 
constants, the response model can be obtained given a stochastic model 
for the input earthquake ground motion. For the case of random 
strength, x, it becomes necessary to determine the distribution function 
F x (x) based on statistical and laboratory tests on individual building 
components as related to the overall structural resistance, e.g., 
Kennedy. 5 Similar tests were used by Freudenthal and Wang to estab- 
lish a representative distribution of the ultimate strength of aircraft 
structures. 6 

In this study, consideration will be limited to the first excursion 
failure only, and the input and response of the structure are both 
treated as random processes. The probability of the duration of the 
response amplitude excursion and other failure mechanisms such as 
wearing and fatigue are not considered. In the first excursion failure, 
a structure is said to have failed if the response parameter, y(t), exceeds 
a prescribed resistance or strength level, x, during the vibration cycles 
caused by the earthquake. Let the duration of the structural vibration 
be approximated by t u , then, for any 1 1 [0, t ] 

p(i) = prob [ | y(t) | ^ x\ ^ t g Q (17) 

- i - wm 

in which W,(t ) = prob [ | y(t) | < x; ^ t ^ t„] = the reliability of 
the system. Two different situations in the structure's resistance char- 
acteristics will be considered below. 

2.4.1 Deterministic Resistance Variable, x 

A random process model for the response of structures subjected to 
a stationary earthquake excitation can be established as follows. A 
simple structure can generally be treated as a lightly damped linear 
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oscillator and its response, y(t), is related by a second-order differential 
equation to the excitation, e.g., a,(0, the ground acceleration of an 
earthquake with intensity i. A multistory structure can be treated 
similarly in generalized coordinates considering normal mode vibrations. 
The structural response y(t) in our case is a Gaussian process which 
approaches stationarity after a few cycles of initial transient motions. 
Let Pi(t)dt denote the probability that y(t) exceeds the threshold y = x 
during the interval [t, t + dt] for the first time since the initial time 
t = 0. The probability density function pi(t), referred to as the first- 
crossing density, is related to the reliability function by — (dW/dt) = 
Pi(t). While establishing the precise behavior of p x (t) for small t poses 
some difficulty, for most practical purposes some approximations can 
be made for large mean failure time t. The simplest approximation to 
the first-crossing density is to assume that the up-crossings of the 
threshold occur rarely in the stationary response, so that they can be 
considered as statistically independent events. If so, the instants at 
which \y(t) | cross the level x from below would constitute a Poisson 
process with a constant rate 2v x , where v T is the level crossing rate of 
y(t) at the level y = x. In this situation it can be easily shown 7 that 

p(i) = 1 - exp (-2v x t ) 

= 1 - exp [-2v t exp (-s 2 /2<r 2 „)] (18) 

in which v„ = ^/(tto-,,) = zero crossing rate of y(t); v x = v u \ v =* = 
v„ exp (—x 2 /2ffl); where a„ = standard deviation of y(t), a* = the 
standard deviation of y(t) = dy(t)/dt. All these quantities are dependent 
on earthquake parameters (therefore, on intensity i) and building 
parameters. These dependences will be derived later in this study. 
More specifically, it will be shown that a y and a 6 are directly proportional 
to i and that i>„ is a constant. 

Expressions for p(i) under other assumptions on the response process 
are presented in Appendix A. 

2.4.2 Random Resistance Variable, x 

The resistance variable x is a random variable with probability 
density function j x {x). In this situation, the failure probability is 
given by 

p(i) = f prob [max \ y(t) | ^ x; ^ t g Qj x {x) dx. (19) 

■'0 t 

Let Y(t B ) = (max, | y(t) \ , ^ t ^ t ); and N„(t) = the number of 
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peaks of | y(t) | in the time I; then \N u (t); g t g t \ is a random 
process. Assuming it is a stationary Poisson process of intensity X„ , 
it follows from eq. (18) that X,. - 2i>, . For situations as described in 
Appendix A, X„ = 2v ln[l - exp {-x 2 /2&j)] from assumption (i); 
and X„ = 2v I \/2£ x/o, l from assumption (ii). 

The first excursion probability can be expressed in terms of X„ as: 

prob [Y(Q £ x] - exp [-X W < (1 - F T (x))} (20) 

in which F x (x) = f_ M j x (x') dx'. From eqs. (19) and (20) the expression 
for p(i) becomes 



v (j) = p ji _ exp [-X U U1 - F x (x))])j x (x) dx. 



(21) 



A closed form solution of eq. (21) is possible when f x (t) has a simple 
expression such as a uniform, Gaussian, or Rayleigh density function. 
In general, eq. (21) can be conveniently solved by numerical integration. 

2.5 Ground Motion Statistics 

The statistics characterizing the random ground motion shall now 
be brought into the formulation. Let a,-(0 be the ground accelerations 
of earthquakes of intensity i and assume {o,-(0; ^ t ^ t \ be a sta- 
tionary process with a power spectral density function G' ,(co), where 
co is the frequency variable. Such a stochastic ground motion model 
has been proposed and used extensively, e.g., by Liu 8 and by Jennings 
et al." Further, assume that the process Ml)) is a filtered white noise 
with a constant power spectrum density G per unit intensity, and that 
the ground filter is a linear, single-mode oscillator with constant fre- 
quency and damping characteristic values u) tt and £„, respectively. The 
following is a derivation of eq. (22) showing the direct relationship 
between (?„,-(«) and intensity i. 

G ai (i0) = Y^ 2x2 T , t 2 2 2' ( 22 ) 

The earthquake intensity value is a measure of the damage potential 
which is represented by the corresponding response spectrum S, = 
max, | y(t) \ for a structure with natural frequency and damping pa- 
rameters Wo and |„ . Therefore, S t = S„(i, u>„ , £„ , t„) is clearly an in- 
creasing function of i. The precise functional relationship between 
S 6 and i is not yet known, but can be obtained from data fittings of 
calculated response spectral values of past earthquakes with known i. 
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A simple linear approximation may be made for our analysis by assuming 

S a (i, « , £„ , = ki (23) 

where k is a constant of proportionality. From the ground acceleration 
model defined above, it can be shown 10 that 

S, = u„Ka u (24) 

where 

K = (2 In v e t ) h + 0.577(2 In M,)"*- (25) 

Since K is independent of i, as will be shown later in this section, 
it is obvious that 

G oi (u>) = i 2 GM (26) 

satisfies eqs. (23) and (24); and according to the ground motion model 
as defined earlier, G„{u>) is given 10 by 

Sofa) = 7"2 2n2 , At 2 2 2' (27) 

Finally, eq. (22) follows directly from eqs. (26) and (27). 

It may be noted from eqs. (23), (24), and the relation o- 2 = /"„ (7 al (a>) 
| H(u>) | 2 du, that the power spectrum density G ai (u>) of the earthquake 
process is proportional to i 2 , which agrees with eq. (26). Housner and 
Jennings 11 have used the relationship G ai (u) = const. *S^ , which also 
leads to our assertion of eq. (26). 

A difficulty exists in determining the value of the constant spectral 
density G corresponding to a unit Modified Mercalli intensity level. 
Because the intensity cannot be precisely related to the earthquake 
waveform parameters such as the amplitude of acceleration, velocity, 
displacement, response spectrum intensity, etc., some normalization 
procedures based on these parameters must be used to determine G. 
For example, a constant power spectral density level for the input 
white noise to the ground filter is determined by matching the cor- 
responding expected velocity spectra of the filter's response to Housner's 
average velocity spectra. 12 

Using the well-known relation a 2 = /"„, G(u) dot, (i.e., the variance 
of a random process is equal to the integral of its power spectral density 
over the entire real line representing frequency), it follows from eq. (24) 
that the variances of a t (t) and d,(<) are respectively o- 2 = i 2 Gir/{2£ g u)l) 
and o- 2 = &) 2 o- 2 . Also, from the relation G y (u) = \ H(u)) | 2 G ai (u) in 
which H(u>) = (to 2 — w 2 — 2#„w co)~" 1 = the transfer function of the 
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simple structure for displacement output y(t) and input a t (t), where 
j represents the complex unit, it can be shown that 1 

ttGB V 



1 l ■ AJLfi - AlAlJ * = *'* (28) 

* GA > V i = e t i (29) 



A,B - A A\ 



in which A = u?ja\ ,A t = 2a) o a>,(£ «, + £„co„), A 2 = a>l + a>* + 4£ £ B w o « g , 
^ 3 = 2(£„a> + £„aO and B = A 3 A 3 - A t . 

From eqs. (28) and (29), the zero crossing rate of y(t) is 



/ 1 f ApAJ f 



(30) 



To show that K is independent of i in eq. (25), it is sufficient to show 
that v„ is likewise independent of i. This is obvious from eq. (30) . 

This completes the discussion on the determination of the failure 
probability of structure p(i). It is shown that eqs. (18), (21), and eqs. 
(37) through (39) in Appendix A define p(i) for various failure mechan- 
isms. Furthermore, substituting eqs. (27) through (30) in the appropriate 
terms in p(i), indicates that p(i) is a function of intensity i. Finally, sub- 
stituting eq. (16) and various expressions for p{i) into eq. (10) deter- 
mines p, which is a function of i„ . 

2.6 The Expected Random Earthquake Loss, e 

The earthquake loss depends on earthquake and building parameters, 
and the extent to which human lives are in danger. The quantity e 
can be determined from statistical data of actual earthquake damage. 
Unfortunately, empirical values of e for different classes of constructions 
are not presently available. It is logical to assume that e is directly 
proportional to the damage potential of earthquakes, therefore, either 
of the following relationships, or their combinations, may be appropriate: 

e = d max o,-(0 (31a) 

e = C 2 T (31b) 

e = C 3 {S a ,orS a , or ft}, (31c) 

in which C, , C 2 , C 3 are the constants of proportionality and S repre- 
sents the expected response spectrum associated with the subscript 
response parameter. Equation (31c) in which e is expressed in terms of 
the expected velocity spectrum ft appears superior to others because 
the effects of the amplitude, duration as well as the frequency char- 
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acteristics of the earthquake accelerogram, are all considered. Thus it 
will be used in the following analysis. Let C 3 = tA f j, where e is the per- 
cent loss in building cost per unit response velocity spectrum and e/ is the 
expected earthquake loss per unit floor area per unit response spectrum. 
It should be noted that no upper bounds for e can be established in situa- 
tions where human lives arc involved. In these situations, the value of e 
would increase for large-occupancy structures such as hospitals, schools 
and office buildings, etc., and decrease for small-occupancy structures 
such as warehouses, machine rooms, unmanned equipment buildings, etc. 
The determination of c, with due considerations to loss of human lives, 
needs actual earthquake life-loss statistics and a mathematical model 
which converts life-loss into dollars. These matters will need further 
studies and more data collection. Clearly, expressing e in terms of 
initial building investment is a convenient way of incorporating all 
possible losses in an earthquake environment. A sensitivity analysis 
for the parameter e should provide some insight to the overall cost 
structure. From eqs. (24) and (31c) 

g-ew&r, ^ (32) 

= eA f jo} Kd u i 
where 

12 

i = X) #'(*)• 
i-i 

It can be noted from the above that e is independent of the design 
intensity i D and the expected service life of a building. This is because 
according to its definition, e is the expected loss associated with a 
"single" random event. The quantity e should not be confused with 
the total expected loss of building (see Appendix B for its derivation) 
throughout the entire service life, which should be expected to increase 
with service life but decrease with design intensity. 

2.7 Earthquake Risk Cost 

The expression for the earthquake loss function K d (i„), for the case 
of deterministic structural resistance and independent crossings of 
response process can now be established. Substituting eqs. (8), (32), and 
(10) [//(*'), p(i) and a v given by eqs. (16), (18) and (28)] into eq. (6) 
leads to: 

K& ) = *hA f ju Ke v g{t)(^) 2 g *-"*■ jjj e-'"" 

•[1 - exp (-2t oVo e' x ' /2e "' i ')) } i = 1, 2, ■•• , 12 (33) 
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in which K and v, are given by eq. (25) and oq. (30) respectively. To show 
cq. (33) is a monotonically decreasing function of i , let 



akA,ja KO u g(t)(^J = Q, , 



£<«p(-g-Qi, 1=5, 2U-M, ^.-r. 



and 



«o 



= z.(t p ). 



(34) 



Equation (33) is rewritten as 
Z£i.) = f^ Zk = T l2 - T im , 

k-i» 

2k = e" Jk [l - exp (-*f r/k ')L 

2 = 0, k, i„ = 1, 2, •■• , 12 

where 7\„ = £k-i z k is a function of z'„ and T 12 = 2", * is a constant. 
Notice in eq. (34) that for 5 and f ^ 0, both exp (-5k) and exp (-f/k 2 ) 
are bounded between zero and unity; furthermore, for m = 0, the 
quantities exp [ — n exp (— f/k 2 )] and 1 — exp [ — /x exp ( — f/k 2 )] are 
also bounded between zero and one. It is apparent that since ^ z k ^ 1 
and 7\ = J£k-i ^ = °> therefore Z u is a monotonically decreasing 
function of i n , as is expected. 

Substituting eqs. (3) and (34) into eq. (1), the final expression for the 
objective function is obtained: 

K{i ) = KXQ + KSo) 

= AXf + <*) + Bt0i*.(O 

in which the first term in the right-hand side increases with i„ and the 
second term decreases with i . The optimum intensity i* is determined 
by setting equal to zero the first derivarive of eq. (35) with respect to 
design intensity ?'„ , i.e., 

dZ„(i„) A,c 



di e fii Q 2 



(36) 



Alternatively, i* can be obtained from eq. (35) by direct, numerical 
evaluation of the function K(i„) for all i„ . 

As an illustration of the presented approach, a hypothetical building 
design problem is numerically analyzed for its optimal design earth- 
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quake intensity. Figure 2 is a plot of i versus the normalized total cost 
K{i )/A f j computed using the following data: The cost data, e = 
1 percent, g(t) = 5 percent; the building data, x = deflection = 0.9 X 
10~ 2 ft; the earthquake data, I = 50 miles, a = 3 earthquakes per 
year, = 2, t = 25 s, t = 40 yrs, c, = 8.16, c 2 = 1.45, c 3 = 2.46, m = 5, 
w„ = 31.4 rad/s, £„ = 0.5, co„ = 3.14 rad/s, £ = 0.5, and d = 40 mi, 
and G »* 1.6 </ 2 s~ 3 . Although these numerical values are used for 
illustration purposes, nevertheless, the earthquake parameters reflect 
realistic data based on seismicity in Southern California. Four different 
seismic protective ratios c/f = 0.01 to 0.04 are considered. The results 
indicate that for this specific design problem, the optimal design in- 
tensity for the building is VII. 
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Fig. 2 — Optimum design earthquake intensity analysis for a sample building. 
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A preliminary sensitivity analysis on this numerical example indicates 
that the convexity of the total cost function becomes more apparent as 
the building resistance parameter x increases (in this case building 
deflection measured in feet). For the problem under consideration, the 
value of x above which the design intensity could be established is 
found to be 0.7 X 10~ 2 feet. For smaller values of x, the cost ratio 
K(i„)/A f f becomes insensitive to design intensity i and ratio c/j as 
the total cost curve becomes flat for i > V. 

III. CONCLUSION 

A simple mathematical approach is presented to determine the 
optimal design intensity of earthquakes for structures. The objective 
function to be optimized is taken as the total construction cost of the 
structure plus the expected cost of earthquake damage throughout 
the entire service life of the structure. For the case of deterministic 
structural resistance and probabilistic structural response with Poisson 
(independent) crossings, the objective function is derived in terms of 
the building and earthquake variables. The optimal design intensity 
can then be determined by minimizing the objective function with 
respect to the intensity variable. Other optimum design variables can 
also be obtained by simply regarding them as the decision variables 
in the objective function and by performing optimization analysis. 
The resulting equations are relatively simple and can be easily handled 
for numerical studies and sensitivity analysis. Generalizations of the 
results for nondeterministic structural resistance and for response 
processes different from those having Poisson crossings are also in- 
dicated. 
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APPENDIX A 

Expressions for p(i) under Different Assumptions 

(i) Independent Peaks — The dispersion in the number of peaks 
of the narrowband response y(t) is neglected and the magnitudes of 
these peaks are assumed to be statistically independent variables, 
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and each having the probability distribution P[\ y | peak < x] = 1 — 
(v z /p ), then 7 

p(i) = 1 - exp [-2v t In (1 - e~" /2 "')]. (37) 

(ii) Independent Envelope Crossings — The crossings of the envelope 
of y(t) are assumed to be independent and in this situation 7 

p(t) = 1 - exp [-2(2fc)*(*/cr f )t. exp {-x*/%fi\ for £ « 1 (38) 

in which £ equals the damping ratio of the structure. 

(in) Two- State Markov-Process Assumption 13 — The successive in- 
tervals that the envelope of y(t) spends above and below the level x 
are assumed to be random variables with exponential distributions. 
In this case 

p(i) = 1 - exp [-(1 - vjv.) exp (-5J,)] (39) 

where 

8, = w x (l - v x / Vo y l 

and n z = /" p(r, r)r dr \ r=x is the envelope crossing rate of y(t). 

APPENDIX B 

Expressions for Failure Probability, pit), and Total Expected Loss of 
Building, D 

Let hit) be the expected loss in case of failure and p(t) be the failure 
probability density function for the building, then 



D 



= f h(t)p(t)dt. (40) 

Jo 



where t D is the service life of the building. From the logic leading to 
eq. (32), and assuming a discount factor of cost g(t) = 1 — (t/t D ), the 
function h(t) can be written as 

h(t) = KtiQiv + 9(t)) (41) 

where rj is percent of construction cost, representing the earthquake 
loss. The failure density p(t) = dP(t)/dt, where P(t) is the failure 
probability of the building and is given by 

Pit) -1-EI prob [N t (t) = n](l - p(i)) n (42) 
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where p(i) is given by eqs. (18) and (21), and N&) equals the total 
(random) number of earthquakes of intensity i in t years. According to 
eq. (7) 

Prob [N t (t) = n] = ^ exp (-AF,<rf)(AF t <rf)" (43) 

where AF, is the probability that given an earthquake occurs, this 
earthquake has an intensity equal to i. For a linear source of earth- 
quakes it follows from eq. (11) that 



AF< = 7 TJ exp(-^ 



1 - exp I -- 



I \ C2' _ \ C2/ _ 



(44) 



Equations (40) through (44) completely define the total loss expectation, 
D, and from above, it is obvious that D increases with t D and decreases 
with i„ . 
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